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The lightest neutralino, the gravitino, and the axino 

Frank Daniel Steffen 1 a 
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Abstract. In supersymmetric extensions of the Standard Model, the lightest neutralino, the grav- 
itino, and the axino can appear as the lightest supersymmetric particle and as such provide a 
compelling explanation of the non-baryonic dark matter in our Universe. For each of these dark 
matter candidates, I review the present status of primordial production mechanisms, cosmological 
constraints, and prospects of experimental identification. 

PACS. 95.35.+d Dark matter - 12.60.Jv Supersymmetric models - 04.65.+e Supergravity 



1 Introduction 

Numerous astrophysical and cosmological considera- 
tions point to the existence of non-baryonic dark mat- 
ter in our Universe [HE]- In fact, based on observa- 
tions of supernovae, galaxy clusters, and the cosmic 
microwave background (CMB), we believe today that 
our Universe is flat with about 76%, 20%, and 4% of 
the critical energy density p c provided in the form of 
dark energy, non-baryonic dark matter, and baryons, 
respectively [3j[4] . A nominal "3c" rangcQ of the dark 
matter density J?dm = Pdm/Pc can be inferred from 
measurements of the CMB anisotropies by the Wilkin- 
son Microwave Anisotropy Probe (WMAP) satellite [3] 



Qi<y h 2 _ n in c+o.02i 
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with h = 0.73±g;g| denoting the Hubble constant in 
units of lOOkmMpc^V 1 . 

Relying on the pieces of evidence, we think that 
a particle physics candidate for dark matter has to 
be electrically neutral, color neutral0 and stable or 
have a lifetime Td m that is not much smaller than the 
age of the Universe today to ~ 14 Gyr. Moreover, the 
species providing the dominant contribution to I?dm 
have to be sufficiently slow to allow for structure for- 
mation. For example, since the neutrinos of the Stan- 
dard Model are too light, J2 t m ^ < °( l eV ) 0, they 
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1 Note that the nominal "3er" range is derived assuming a 
restrictive six-parameter "vanilla" model. A larger range is 
possible — even with additional data from other cosmologi- 
cal probes — if the fit is performed in the context of a more 
general model that includes other physically motivated pa- 
rameters such as a nonzero neutrino mass [5]. Thereby, the 
ranee 0.094 < O dln h 2 < 0.136 has been obtained in Ref. [5]- 

2 A colored dark matter candidate is disfavored by severe 
limits from searches for anomalous heavy nuclei [4]. 



were too fast at early times. Accordingly, they are clas- 
sified as hot dark matter which can constitute only a 
minor fraction of J7dm since otherwise structure forma- 
tion cannot be understood [7j. Thus, the observation- 
ally inferred dark matter density can be considered as 
evidence for physics beyond the Standard Model. 

Supersymmetric (SUSY) extensions of the Stan- 
dard Model are an appealing concept because of their 
remarkable properties, for example, with respect to 
gauge coupling unification, the hierarchy problem, and 
the embedding of gravity [8llMT0l[iTl[l2lfT%] . As super- 
partners of the Standard Model particles, new parti- 
cles appear including fields that are electrically neutral 
and color neutral. Since they have not been detected 
at particle accelerators, these sparticles must be heavy 
or extremely weakly interacting. 

Because of the non-observation of reactions that vi- 
olate lepton number L or baryon number B, it is often 
assumed — as also in this review — that SUSY theories 
respect the multiplicative quantum number 



R=(-l) 



3B+L+2S 



(2) 



known as R-parity, with S denoting the spin. Since 
Standard Model particles and superpartners carry re- 
spectively even (+1) and odd (-1) R-parity, its con- 
servation implies that superpartners can only be pro- 
duced or annihilated in pairs and that the lightest su- 
persymmetric particle (LSP) cannot decay even if it is 
heavier than most (or all) of the Standard Model parti- 
clesH An electrically neutral and color neutral LSP can 
thus be a compelling dark matter candidate. For the 
lightest neutralino, the gravitino, and the axino, which 
are well-motivated LSP candidates, this is shown be- 
low. For each scenario, I will address implications for 
cosmology and experimental prospects. Note that the 



While R-parity conservation is assumed in this review, 
its violation is a realistic option; see, e.g., [14 15 , 16 , 171118]. 
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discussion of gravitino/axino dark matter in Sects. [3] 
and[?]will be more extensive than the one of neutralino 
dark matter in Sect. [2] for which numerous excellent 
reviews exist such as |19 p i2 | l20 l l2T ] . 

2 Neutralino Dark Matter 

The lightest neutralino Xi appears already in the min- 
imal supersymmetric Standard Model (MSSM) as the 
lightest mass eigenstate among the four neutralinos be- 
ing mixtures of the bino B, the wino W, and the neu- 
tral higgsinos H® and H®. Accordingly, \i is a spin 1/2 
fermion with weak interactions only. Its mass to^o de- 
pends on the gaugino mass parameters Mi and M<x, on 
the ratio of the two MSSM Higgs doublet vacuum ex- 
pectation values tan/3, and the higgsino mass param- 
eter /i. Expecting m^o = 0(100 GeV), Xi is classified 
as a weakly interacting massive particle (WIMP). 

Motivated by theories of grand unification and su- 
pergravity [22j and by experimental constraints on fla- 
vor mixing and CP violation [4], one often assumes 
universal soft SUSY breaking parameters at the scale 
of grand unification M GUT ; cf. [TTlfTSTfrgpU] and ref- 
erences therein. For example, in the framework of the 
constrained MSSM (CMSSM), the gaugino masses, the 
scalar masses, and the trilinear scalar interactions are 
assumed to take on the respective universal values 
mi/2, m , and A at Mgut- Specifying m 1 / 2 , m , A , 
tan/3, and the sign of /i, the low-energy mass spec- 
trum is given by the renormalization group running 
from Mqut downwards. 

Assuming Aq — for simplicity, the lightest Stan- 
dard Model superpartner — or lightest ordinary super- 
partner (LOSP) — is either the lightest neutralino Xi or 
the lighter stau t\ , whose mass is denoted by . If 
the LSP is assumed to be the LOSP, the parameter re- 
gion in which < m^o is usually not considered be- 
cause of severe upper limits on the abundance of stable 
charged particles [4]. However, in gravitino/axino LSP 
scenarios, in which the LOSP is the next-to-lightest 
supersymmetric particle (NLSP), the Ti LOSP case 
is viable and particularly promising for collider phe- 
nomenology as will be discussed in Sects. [3] and |H 

In Fig. [T] (from j23j) the dotted (blue in the web ver- 
sion) lines show contours of tolosp m the (mi^jWo) 
plane for Aq = 0, /i > 0, tan/3 = 10. Above (be- 
low) the dashed line, m^o < {m^ < m^o). The 
medium gray and the light gray regions at small m 1 / 2 
are excluded respectively by the mass bounds rn~± > 
94 GeV and tojj > 114.4 GeV from chargino and 
Higgs searches at LEP [3]. It can be seen that m^o — 
O(100 GeV) appears naturally within the CMSSM. 

2.1 Primordial Origin 

The Xi' s were in thermal equilibrium for primordial 
temperatures of T > Tf ~ m^o/20. At Tf, the an- 
nihilation rate of the (by then) non-relativistic Xi s 
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Fig. 1. Contours of ttilosp (dotted blue lines) and 1los P 
(solid black lines) in the (mi/2,mo) plane for Ao = 0, 
/i > 0, tan/3 = 10. Above (below) the dashed line, 
m~o < rrir 1 (nif 1 < rn~o). The medium gray and the light 
gray regions show the LEP bounds m~± > 94 GeV and 

mu > 114.4 GeV, respectively [4|. The contours are ob- 
tained with the spectrum generator SuSpect 2 . 34 [24] us- 
ing m t = 172.5 GeV and m b (m b ) MS ' = 4.23 GeV, and with 
micrOMEGAs 1.37 [25lf26] ■ From [23]. 

becomes smaller than the Hubble rate so that they 
decouple from the thermal plasma. Thus, for T < Tf, 
their yield Y^o = n^o/s is given by Yij, 00 » Y~f(T f ), 

where n~ q ' is the (equilibrium) number density of xV s 

and s = 2ir 2 g*s T 3 /45 the entropy density. Depend- 
ing on details of the Xi decoupling, Y"~ OC is very sen- 

sitive to the mass spectrum and the couplings of the 
superparticles. Indeed, convenient computer programs 
such as DarkSUSY [S7J or micrOMEGAs 1 . 37 [2%]l2l)] are 
available which allow for a numerical calculation of the 
LOSP decoupling and the resulting thermal relic abun- 
dance in a given SUSY model. 

The F^osp contours shown by the solid black lines 
in Fig. Q] illustrate that the Xi LSP yield can easily 
vary by more than an order of magnitude. Because of 
this sensitivity, the associated thermal relic density 

Q^h 2 = m^o Y^ c s(T ) h 2 / Pc (3) 

agrees with Q^Ji only in narrow regions in the pa- 
rameter space; p c /[s(T Q )h 2 } = 3.6 x 10" 9 GeV 0]. This 
can be seen in Fig. [2] (from [28 ) where the black strips 
indicate the region with 0.087 < %oft 2 < 0.138. 

Remarkably, it is exactly the small width of the 
i?£0 = f?dm regions which could help us to identify 
Xi dark matter. Once sparticles are produced at col- 
liders, the data analysis will aim at determinig the 
SUSY model realized in nature [2"9"1[3"U] . For the recon- 
structed model, a precise calculation of fl^o is possible 
assuming a standard thermal history of the Universe. 
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Fig. 2. Regions (black) with 0.087 < %o/i 2 < 0.138 in 
the (mi/2, m o) plane for Aq = 0, fj, > 0, tan/3 = 10, 
and rat = YJ2.7 GeV. In the dark gray triangular re- 
gion, m~o > m,r x . The light gray region at small mi/2 
is excluded by the requirement of correct electroweak sym- 
metry breaking or by sparticle search limits [28], the two 
medium shaded (light pink in the web version) bands by 
the LEP bound mu > 114 GeV, and the small light shaded 
(green in the web version) spot by the b — > sj constraint: 
2.65 < BR(6 -> s 7 )/10" 4 < 4.45. The dark shaded (red in 
the web version) band is compatible with having a Stan- 
dard Model like Higgs boson near 115 GeV. From [28] . 



Because of the sensitivity of f2^o with respect to the 
SUSY model, an agreement of the obtained fi^o with 
/2dm will then be strong evidence for the X? LSP pro- 
viding fidm and for a standard thermal history up to 
the x?-decoupling temperature Tf . Since xV s decouple 
already as a non-relativistic species, it is also guaran- 
teed that they are sufficiently cold to allow for cosmic 
structure formation. 



2.2 Experimental Prospects 

For experimental tests of the X? dark matter hypothe- 
sis, three complementary techniques exist: indirect, di- 
rect, and collider searches. While there is an enormous 
activity in each of those fields, I will summarize only 
the main ideas. For more detailed discussions, see [311 
|2"TII3"2] and references therein. 

Let us first turn to indirect searches. Since dark 
matter clumps, one expects regions with an increased 
Xi density such as galaxy halos, the center of galax- 
ies, and the center of stars. While Xi P&i r annihilation 
after Xi decoupling is basically negligible for calcu- 
lations of fi^a , it should occur at a significant rate 
in these regions. The resulting Standard Model par- 
ticles should then lead to energetic cosmic rays and 
thereby to an excess of photons, neutrinos, positrons, 



and antiprotons over backgrounds expected from stan- 
dard cosmic ray models without dark matter annihi- 
lation. In fact, data from the Energetic Gamma Ray 
Experiment Telescope (EGRET) has already been in- 
terpreted as evidence for Xi annihilation [33[134] within 
SUSY models that will be testable in direct and col- 
lider searches. For a discussion of these and other po- 
tential hints, see [5SlT2Tj and references therein. 



In direct searches, one looks for signals of xV s — or 
more generally WIMPs — passing through earth that 
scatter elastically off nuclei. Being located in environ- 
ments deep underground that are well shielded against 
unwanted background, an enormous sensitivity has al- 
ready been reached [36 ,37 ,38, 39]. Since no unambigu- 
ous signal of a Xi _nuc l eus scattering event has been 
observed so far, m^o -dependent upper limits on the 
respective Xi cross section are obtained. Indeed, the 
current best limits given by the CDMS II [38 and 
the Xenon 10 |39j experiments exclude already a part 
of the SUSY parameter space; see, for example, [2"0l 
[2T1I32] and references therein. These limits, however, 
depend on the assumed Xi nux at the detector loca- 
tion, which is subject to significant uncertainties due 
to possible inhomogeneities in the dark matter distri- 
bution in galaxies. Such inhomogeneities should mani- 
fest themselves also in indirect searches which can help 
to reduce those uncertainties. Once Xi events are ob- 
served in direct searches, one can succeed in recon- 
structing the Xi velocity distribution [40] . By analyz- 
ing the recoil spectra, m^o can even be estimated in a 
way that is independent of the dark matter density on 
earth HTl. 



In most searches for SUSY at colliders, it is as- 
sumed that R-parity is conserved. Accordingly, one 
expects that superpartners are produced in pairs be- 
fore decaying via cascades into the LSP and energetic 
fermions. As a weakly-interacting particle, every X? 
LSP produced will escape the detector without leav- 
ing a track. Thus, the existence of SUSY and the Xi 
LSP has to be inferred from studies of missing trans- 
verse energy E^ iss and energetic jets and leptons emit- 
ted along the cascades. Along these lines, ongoing in- 
vestigations are pursued based on data from pp col- 
lisions with a center-of-mass energy of s/s = 2 TeV 
at the Fermilab Tevatron Collider. While lower lim- 
its on the masses of squarks and gluinos have been 
extracted, no evidence for SUSY or the X? LSP has 
been reported so far |42,43J. With the first pp colli- 
sions with y/s — 14 TeV at the CERN Large Hadron 
Collider (LHC) expected in the year 2008, there are 
high hopes that the new energy range will allow for a 
copious production of superpartners. Here large E™ 1BS 
will be the key quantity for early SUSY searches [H] 
145] . Despite an enormous potential for mass and spin 
measurements of SUSY particles at the LHC |46j , ad- 
ditional precision studies at the planned International 
Linear Collider (ILC) [47[l48] appear to be crucial for 
the identification of the x? LSP [3TU49] . 
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3 Gravitino Dark Matter 

The gravitino G appears (as the spin-3/2 superpartner 
of the graviton) once SUSY is promoted from a global 
to a local symmetry leading to supergravity [8j. The 
gravitino mass depends strongly on the SUSY- 
breaking scheme and can range from the eV scale to 
scales beyond the TeV region [9l[TT|[50l[5Tl[52p3l[54l 
[55] . For example, in gauge-mediated SUSY breaking 
schemes 50. 5"Tll5*2] . the mass of the gravitino is typi- 
cally less than 1 GeV, while in gravity-mediated sche- 
mes [9]ITT] it is expected to be in the GeV to TeV range. 
The gravitino is a singlet with respect to the gauge 
groups of the Standard Model. Its interactions — given 
by the supergravity Lagrangian 56,5] — are suppressed 
by the (reduced) Planck scale [J| 



M P = 2.4 x 10 18 GeV. 



(4) 



Once SUSY is broken, the extremely weak gravitino in- 
teractions are enhanced through the super-Higgs mech- 
anism, in particular, at energy/mass scales that are 
large with respect to mg. Nevertheless, the gravitino 
can be classified as an extremely weakly interacting 
particle (EWIP). It must not be massive since even 
a light gravitino can evade its production at colliders 
because of its tiny interaction strength. Considering 
the case of the G LSP, in which the LOSP is the un- 
stable NLSP that decays eventually into the G LSP, 
< m^o (cf. Fig.QJ is viable as already mentioned^ 



3.1 Primordial Origin 

Assuming that inflation governed the earliest moments 
of the Universe, any initial population of gravitinos 
must be diluted away by the exponential expansion 
during the slow-roll phase. Indeed, gravitinos are typ- 
ically not in thermal equilibrium with the primordial 
plasma after inflation because of their extremely weak 
interactions^ At high temperatures, however, they can 
be produced efficiently in thermal scattering of parti- 
cles in the primordial plasma. Derived in a consistent 
gauge-invariant treatment, the resulting thermally pro- 
duced (TP) gravitino density reads [64]l6"5] 




(5) 



gi are understood to be evaluated at the reheating tem- 
peraturejj 1 after inflation Tr [55] For the case of uni- 
versal Mj.,2,3 = m-i/2 at Mqut and mg <C Mi, i.e., 
(l + M?/3m| 

mated by the convenient expression [67 



M 2 /3to~, L?~ p /i 2 can be approxi- 

' Cx Cr 



0.32< 



/ 10 GeV \ / toi/ 2 y / i R \ 
V m d )\1 TeV/ VlO 8 GeV/' W 



The thermally produced gravitinos do not affect the 
thermal evolution of the LOSP (or NLSP) prior to 
its decay which occurs typically after decoupling from 
the thermal plasma. Moreover, since each NLSP de- 
cays into one G LSP, the NLSP decay leads to a non- 
thermally produced (NTP) gravitino density [68.69, 

mm 



f2f p h 2 = m d Y^ SpS (T a )h 2 /p, 
so that the guaranteed density is given 



n s h 2 = Hfh 2 



a 



NT P/l 2 

Cr 



(7) 



(8) 



While /2~ p is sensitive to Mi and Tr for a given 

Cr 

mg, /2~ TP depends on Ynu3 P = ^losp anc ^ thereby 
on details of the SUSY model realized in nature; cf. 
Sect. 12.11 For the case of the T\ NLSP, simple approx- 
imations can be used such as [69ll 7411751 



Y- 



dec 



0.7 x 10 



-12 



1 TeV 



(9) 



which is valid outside of the T\-Xi coannihilation re- 
gion for a spectrum in which is significantly be- 
low the masses of the lighter selectron and the lighter 

~ B with a 
- i^dni are 



smuon, m?! 
mass of m 



*C TO ei,pii and in which \i 
= 1.1 rtif x LJ Scenarios with Qq 



6 For a discussion on the Xr definition, see Sec. 2 in [23] . 

7 Note that the field-theoretical methods applied in the 



derivation of ([5} 64 65 require weak couplings ji<l and 
thus T ^> 10 6 GeV. For an alternative approach, see [66] . 

8 In this review I do not discuss gravitino production 
from inflaton decays which can be substantial depending 
on the inflation model; see, e.g., |72II73| . 

9 The Flosp contours in the ri LOSP region in Fig. [1] in 
which rriT X J; fn^ 1 ^ 1 < 1.1 m^, illustrate that the t\ LSP 
yield can be about twice as large for a given due to 
slepton coannihilation. Approaching the Xi-ri coannihila- 
tion region, m~o « , even larger factors occur. 



Table 1. Assignments of the index i, the gauge coupling gi, 
and the gaugino mass parameter Mi to the gauge groups 
U(1)y, SU(2)l, and SU(3) C , and the constants u>i and hi. 



with u>i, the gauge couplings gi, the gaugino mass pa- 
rameters Mi , and ki as given in Table Q] Here Mi and 



4 A stop h NLSP is not feasible in the CMSSM [57]. 

5 In gauge-mediated SUSY breaking scenarios, light 
gravitinos can be viable thermal relics if their abundance is 
diluted by entropy production, which can result, for exam- 
ple, from decays of messenger fields 58 , 59 , 60 , 61 , 62 , 63 . 
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Fig. 3. Upper limits on the reheating temperature Tr in 
the G LSP case. On the upper (lower) gray band, fl-^ G 
Qll, for Mj.,2,3 = mi/ a = 500 GeV (2 TeV) at M GUT . The 
corresponding limits from i?~ p /i 2 < 0.126 shown by the 
dashed and dotted lines are obtained respectively with ((5} 
for Mi/10 = M 2 /2 = Mz = mi/a at Mgut and with the 
result of Ref. [64] for M 3 = m 1/2 at M GU t- From [23]. 



found for natural mass spectra and for a wide range of 
ttIq—Tr combinations. This is illustrated in Figs. [3] [4] 
and El 



In G LSP scenarios, upper limits on Tr can be 
derived since H% p < Sl dm [7511551177117511731125] . These 
mg-dependent limits are shown in Fig. [3J (from [23]) 
and can be confronted with inflation models. More- 
over, Tr limits are important for our understanding 
of the baryon asymmetry and, in particular, for ther- 
mal leptogenesis 179 1 180 ) . For given f2% p , the bound 



, Sl~ gives upper limits on and • 



The limits obtained for a t% NLSP with ([9]) are shown 
in Fig. SI In Fig. El (from [8j) regions with % <E 
are shown for T R = 10 7 , 10 8 , and 10 9 GeV. Here both 
f?~ p and ,T2~ TP are taken into account for = mo 
within the framework of the CMSSM. 



While thermally produced gravitinos have a neg- 
ligible free-streaming velocity today, gravitinos from 
NLSP decays can be warm/hot dark matter. In the 
T\ NLSP case, for example, upper limits on the free- 
streaming velocity from simulations and observations 
of cosmic structures exclude < 0.7 TeV for i?~ TP ~ 
■^dm [7J] • Such scenarios (gray band in Fig. UJ, how- 
ever, require > 0.7 TeV anyhow and could even 
resolve the small scale structure problems inherent to 
cold dark matter [811871188] . 




10 20 



50 100 200 500 1000 
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Fig. 4. Cosmological constraints on the masses of the grav- 
itino LSP and a purely 'right-handed' ri NLSP. The gray 



3<x 
dm ■ 



Above this band, ftp > 0.126. 



band indicates J?~ TP £ SI, 
Only 10% (1%) of ft m is provided by J?g TF for scenar- 
ios that fall onto the thin solid line labeled by / = 0.1 
(0.01). The dotted lines show contours of t^ x . The region 
below the long-dash-dotted (red in the web version) line 
and below the long-dashed (green in the web version) line 
is disfavored by the observationally inferred abundances of 
primordial 6 Li [67] and 3 He/D [ST]. The effect of electro- 
magnetic and hadronic energy injection on primordial D 
disfavors the regions inside the short-dash-dotted (blue in 
the web version) curves and to the right or inside of the 
short-dashed (blue in the web version) curves, respectively. 
With ^ and e cm = 0.3E T , the curves are obtained from 
the severe and conservative upper limits defined in Sec. 4.1 
of [74] based on results from [551151] . 



3.2 Cosmological Constraints 

In the G LSP case with conserved R-parity, the NLSP 
can have a long lifetime tnlsp E3 This is illustrated by 
the dotted tnlsp contours in Figs. [JJ and El In partic- 
ular, for the Ti NLSP, one finds in the limit m T — > 0, 



r- 1 {r 1 ^Gr) = 



48ttto 2 ~M| 



m| 



(10) 

while the expression for the lifetime of the Xi NLSP is 
given in Sec. IIC of Ref. [71]. 

If the NLSP decays into the G LSP occur during 
or after big-bang nucleosynthesis (BBN) , the Standard 
Model particles emitted in addition to the gravitino 
can affect the abundances of the primordial light ele- 
ments. Indeed, these BBN constraints disfavor the Xi 
NLSP for m d > 100 MeV [7lT|77||75||55] . For the slep- 
ton NLSP case, the BBN constraints associated with 
electromagnetic/hadronic energy injection have also 
been considered and found to be much weaker but still 



For the case of broken R-parity, see, e.g., [16U17lfT8] . 
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Fig. 5. CMSSM regions with Qgh 2 G fl^ for T R = 10 7 , 
10 s , and 10 9 GeV indicated respectively by the light, 
medium, and dark shaded (green in the web version) bands 
in the (m 1 / 2 ,mo) planes for tan/3 = 10, Ao — 0, /i > 0, and 
mg = mo. The regions excluded by the chargino and Higgs 
mass bounds and the line indicating rn~o = are identi- 
cal to the ones shown in Fig. [T] In the dark gray region, the 
gravitino is not the LSP. The dotted lines show contours of 
the NLSP lifetime. The region to the left of the long-dash- 
dotted (red in the web version) line and to the left of the 
thin gray (pink in the web version) line is disfavored by the 
observationally inferred abundances of primordial 6 Li [67] 
and 3 He/D [81]. The effect of hadronic energy injection on 
primordial D |74] disfavors the t\ NLSP region above the 
short-dash-dotted (blue in the web version) lines. The x? 
NLSP region is disfavored by BBN constraints from energy 
injection [84ll71H77ll78ll85] . On the solid vertical line (violet 
in the web version) nig — 2.5 TeV. From [83] . 



significant in much of the parameter space [71 ] l77 | l78 l 
[71] as can be seen in Fig. [31 where the constraints from 
electromagnetic and hadronic energy release are shown 
respectively by the short-dashed and the short-dash- 
dotted (blue in the web version) lines. The hadronic 
constraints are also shown in Fig. [5] 

It has been realized only recently that already the 
mere presence of long-lived negatively charged parti- 
cles can affect BBN substantially via bound-state ef- 
fects PPPTP5 1ig21ip|IMP 5 p g] . In particular, bou- 
nd-state formation of r{~ with 4 He can lead to an over- 
production of 6 Li via the catalyzed BBN (CBBN) re- 
action ( 4 HeA~) + D ->■ 6 Li + X~ 89J. Thereby, the 
observationally inferred upper limit on the primordial 
6 Li abundance 6 Li/H| obs < 2 x KT 11 [82], leads to the 
CBBN constraint shown by the long-dash-dotted (red 
in the web version) lines in Figs. [4] and (5[ as obtained 
from [67]. Indeed, for a typical yield ([9]), the 6 Li/H| b s 
limit quoted above implies the constraint [Mll92U9.il 
E7J67]: t ?1 < 5 x 10 3 s. While numerous other CBBN 
reactions can affect the abundances of 6 Li and other 
primordial elements significantly [85 ,93,96, 98] , the ap- 



proximate t^j bound is relatively robust. In fact, by 
systematically taking into account the uncertainties in 
the relevant nuclear reaction rates, it is shown explic- 
itly in Fig. 14 of [96] and in Fig. 5 of [98] that cosmo- 
logically allowed regions for > 10 5 s are extremely 
unlikely. In particular, the 3 He/D constraint on elec- 
tromagnetic energy release 99 becomes severe and can 
exclude t ?1 > 10 6 s 78,85,94,96 . This is shown by the 
long-dashed (green in the web version) line in Fig. [4] 
and by the thin gray (pink in the web version) line in 
Fig. [5] which are obtained from Fig. 42 of Ref. [81] for a 
'visible' electromagnetic energy of E vis = e em = 0.3 E T 
of the tau energy E T = (m? i -m|+mj)/2mr 1 released 

in ri -> GrE3 

The observed Planck spectrum of the cosmic mi- 
crowave background (CMB) provides an additional con- 
straint |100pi01j which is not shown. Indeed, the CMB 
limit derived in |101j is everywhere less severe than the 
severe electromagnetic limit given by the short- 

dashed (blue in the web version) line in Fig. [4] 

As can be seen in Fig. dj the cosmological con- 
straints provide an upper bound on mg once mf t is 
measured. This bound implies upper bounds on the 
SUSY breaking scale, L?g TP , and T R . 

Figure [5] shows that the cosmological constraints 
imply a lower limit on mi/2 |85[[23i and an upper limit 
on Tr [23]. Indeed, from t? x < 5 x 10 3 s, mg-dependent 
limits on the gaugino mass parameter, 

/ \2/5 
mi / 2 >°-9Tev(7o^y . (ID 

and the reheating temperature, 

T R <4.9xl0 7 Gev( I ^| v ) 1/5 , (12) 

have been derived within the CMSSM [STJ E3 While the 
Tr bound can be restrictive for models of inflation and 
baryogenesis, the bound can have implications 

for SUSY searches at the LHC. Depending on mg, 
([TT]) implies sparticle masses which can be associated 
with a mass range that will be difficult to probe at the 
LHC. This is illustrated by the vertical (violet in the 
web version) line in Fig. [5] which indicates the gluino 
mass m-g = 2.5 TeV [83] PI 

3.3 Experimental Prospects 

Because of its extremely weak couplings, gravitino dark 
matter is inaccessible to direct and indirect searches if 
R-parity is conserved0 Also the direct production of 

With a finely tuned m^-mg degeneracy leading to 
-E v i s — > can any bound on energy release be evaded. 

12 Similar limits have recently been discussed in models 
where the ratio mg/mi« is bounded from below [102] . 

13 Note that the mass of the lighter stop is mq 1 — O.lm-g 
in the considered ri NLSP region with mh > 114.4 GeV. 

14 For broken R-parity, G dark matter is unstable so that 
decay products can appear in indirect searches [1 TU 1031 
[104] . 
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gravitinos at colliders is strongly suppressed. Instead, 
one expects a large sample of (quasi-) stable NLSPs if 
the NLSP belongs to the MSSM spectrum. 

In the Ti NLSP case, each heavier superpartner 
produced will cascade down to the T\ which will ap- 
pear as a (quasi-) stable particle in the detector. Such 
a heavy charged particle would penetrate the collider 
detector in a way similar to muons [105, 106, 107J. If the 
produced staus are slow, the associated highly ionizing 
tracks and time-of-flight measurements will allow one 
to distinguish the t\ from a muon 105, 106, 107, 108 . 
With measurements of the t\ velocity (3^ = / c and 
the slepton momentum = \Pt x \, m?! can be de- 
termined: m fl = pfi(l - /32) 1/2 //3fi [105] . For the 
upcoming LHC experiments, studies of hypothetical 
scenarios with long-lived charged particles are actively 
pursued [109 l lllO p ill | flT2]. For example, it has been 
found that one should be able to measure the mass 
of a (quasi-) stable t\ quite accurately [109|J110j P^l 

If some of the staus decay already in the collider de- 
tectors, the statistical method proposed in [108] could 
allow one to measure the T\ lifetime. With (flO]) and 
the measured value of , one will then be able to 
determine also the gravitino mass mg and thereby the 
scale of SUSY breaking. As a test of our understanding 
of the early Universe, it will also be interesting to con- 
front the experimentally determined (rag , ) point 
with the cosmological constraints in Fig. [4J 

Ways to stop and collect charged long-lived parti- 
cles for an analysis of their decays have also been pro- 
posed [H4l[n^[n6l[lTn[n8l[ll9] . It was found that up 
to O(10 3 -10 4 ) and O(10 3 -10 5 ) n's can be trapped per 
year at the LHC and the ILC, respectively, by placing 
1-10 kt of massive additional material around planned 
collider detectors |115pil6| . A measurement of can 
then be used to determine as already described 
above. If mg can be determined independently from 

the kinematics of the 2-body decay t\ Gt, 



m G ~ V m ?i + m T ~ 2m fl E T , (13) 



(see Fig. |4|03 Accordingly, alternative methods such 
as the ones proposed in [12411125) could become essen- 
tial to identify the gravitino as the LSP. 



4 Axino Dark Matter 

The axino a [T2^1ll27U1281fT2^] appears (as the spin- 
1/2 superpartner of the axion) once the MSSM is ex- 
tended with the Peccei-Quinn mechanism 130, 131] in 
order to solve the strong CP problem. Depending on 
the model and the SUSY breaking scheme, the ax- 
ino mass rria can range between the eV and the GeV 
scale |128lll32l[T3"31ll341ll35ill36j . The axino is a sin- 
glet with respect to the gauge groups of the Standard 
Model. It interacts extremely weakly since its cou- 
plings are suppressed by the Peccei-Quinn scale [4] 
H37lll38U139j 

f a > 5 x 10 9 GeV (15) 

and thus can be classified as an EWIP. The detailed 
form of the axino interactions depends on the axion 
model under consideration. We focus on hadronic (or 
KSVZ) axion models p^OllMl] in a SUSY setting, in 
which the axino couples to the MSSM particles only 
indirectly through loops of additional heavy KSVZ 
(s)quarks. Considering a LSP scenarios in which the 
LOSP is the NLSP, < m^o is again viable. 

Before proceeding, it should be stressed that the 
bosonic partners of the axino, the axion and the sax- 
ion, can have important implications for cosmology: 
(i) The relic density of axions can contribute signifi- 
cantly to the dark matter density 137,138,142] and 
thereby tighten the constraints from fi~ p < f2^ m dis- 
cussed below, (ii) Late decays of the saxion can lead to 
significant entropy production [143[|144[|145l H46 and 
thereby affect the cosmological constraints |147j . In 
this review, however, a standard thermal history is as- 
sumed which implies that saxion effects are negligible. 



the lifetime 



can allow for a measurement of the 
Planck scale pplll2lllll8lll22j 



M| 



487T m~ 



m 2 fi 



(14) 



An agreement with which is inferred from New- 
ton's constant [4] G N = 6.709 x 10" 39 GeV~ 2 , would 
then provide strong evidence for the existence of super- 
gravity in nature 120 . In fact, this agreement would 
be a striking signature of the gravitino LSP. Unfortu- 
nately, the required kinematical determination of mg 
appears to be feasible only for [121[|118[|122] mg /uif 1 > 
0.1 which seems to be disfavored according to our 
present understanding of the cosmological constraints 



15 (Quasi-) stable ti's could also be pair-produced in in- 
teractions of cosmic neutrinos in the earth matter and be 
detected in a neutrino telescope such as IceCube [113] . 



4.1 Primordial Origin 

Because of their extremely weak interactions, the tem- 
perature Tf at which axinos decouple from the thermal 
plasma in the early Universe is very high. For exam- 
ple, an axino decoupling temperature of Tf « 10 9 GeV 
is obtained for f a = 10 11 GeV [l33l[l48] . For T R > T f , 
axinos were in thermal equilibrium before decoupling 
as a relativistic species so that [133,149,150,148 



2keV 



(16) 



For Tr < Tf, axinos are not in thermal equilibrium 
with the primordial plasma after inflation but can be 



16 Note that the cosmological constraints described in 
Sect. 13.21 assume a standard thermal history. In fact, en- 
tropy production after NLSP decoupling and before BBN 
can weaken the BBN constraints significantly [1231123] . 



Frank Daniel Steffen 



Review 



generated efficiently in scattering processes of parti- 
cles that are in thermal equilibrium with in the hot 
SUSY plasma [M9l[T50lfT48] . Within SUSY QCD, the 
associated thermally produced (TP) axino density can 
be calculated in a consistent gauge-invariant treatment 
that requires weak couplings (g s <C 1) |148j : 



r „ ,'1.108\ / 10^GeV \ 2 
Tr 



0.1 GcV J \ 10 4 GeV 



(17) 



with the axion-model-dependent color anomaly N of 
the Peccei-Quinn symmetry and the strong coupling 
g s understood to be evaluated at Tr. The thermally 
produced axinos do not affect the thermal evolution of 
the LOSP (or NLSP) which decays after its decoupling 
into the a LSP. Taking into account the non-thermally 
produced (NTP) density from NLSP decays [T5T1H50] 



NTP i 2 



ydec „^rp ^2 , 



l a * NLSP 



s(T Q )h /p c , 



the guaranteed axino density if 



;i7i 



(18) 



(19) 



In Fig. © (from 148) the (m 5 ,T R ) region with 
0.097 < Qj v < 0.129 for f a /N = 10 11 GeV is shown 
by the gray band. Note that (jTTJ) shows a different 
dependence on the LSP mass than the correspond- 
ing expression in the G LSP case (0). Accordingly, 
one finds the different tolsp dependence of the Tr 
limits inferred from &~Jq < fldm- Since thermally 

produced axinos are generated in kinetic equilibrium 
with the primordial plasma, they have a thermal spec- 
trum which allows for the m^-dependent classification 
into cold, warm, and hot dark matter [150j shown in 
Fig. O As can be seen, the Tr limit does not exist for 
mjr < 0.2 keV because of the equality of a production 
and a disappearance rates for T > Tf w 10 9 GeV. With 
a thermal relic density (|16|) in this regime, there will 
be a limit on to^ depending on the constraints inferred 
from studies of warm/hot dark matter [7]. 

The non-thermally produced axino density Q~ TP 
differs from the corresponding expression in the G LSP 
case ([7]) only by the obvious difference in tti~,q. In 

particular, for given Q~ F ', the bound ft~ TP < Odm — 

f2~ q ^ TP as obtained with ([9]) implies limits on to^ and 
which can be read off directly from Fig.|4]after the 
replacement — > nia- Note, however, that the 
contours and the cosmological constraints arc different 
in the axino LSP case. For the T\ NLSP, the following 
lifetime was estimated [124] 



r ?1 ~ T" 1 ^! -► to) ~ 25 s C 



1 - 



( 100 GeV\ / / Q /C a YY \Y 100 GeV x2 



10 1 



GeV/ 



. (20) 



10 



10' 



> 10 
u 




Eh 



10' 



10" 



10 



j I I I llllj 1 I I I llllj 1 — I I 1 1 1 III 1 I [I 111 



I I 



f a /N = 10" GeV- 



OWMAPl2 _ n iiq+0.016 
CDM " — u ' llc> -0.018 



excluded 



hot 



warm 
I i 



cold 
I I 



10 10 10"" 

m a [GeV] 

Fig. 6. Upper limits on the reheating temperature Tr in 
the a LSP case for f a /N = 10 11 GeV. On (above) the gray 
band, Qj v h 2 G 0.113±o.oi| (n~ P h 2 > 0.129). Thermally 
produced axinos can be classified as hot, warm, and cold 
dark matter 150 as indicated. From 148 . 



Axino production in inflaton decays is not considered. 



where the KSVZ-model dependence is expressed by 
C a YY — 0(1) and the uncertainty of the estimate is 
absorbed into £ ~ 0(1). One thus finds a ri lifetime 
in the a LSP case that cannot be as large as the one in 
the G LSP case (p~0|) . Accordingly, the BBN constraints 
are much weaker for the a LSP. For discussions of a 
LSP scenarios with the x? NLSP, see [T5Tlll50lll52j . 
For both the t x NLSP and the Xi NLSP, it has been 
shown that non-thermally produced axinos with mj < 
10 GeV would be warm/hot dark matter [88] . 



4.2 Experimental Prospects 

Being an EWIP, the axino LSP is inaccessible to any 
direct and indirect dark matter searches if R-parity 
is conserved. Also the direct a production at collid- 
ers is strongly suppressed. Nevertheless, (quasi-) stable 
Ti's could appear in collider detectors (and neutrino 
telescopes |113j ) as a possible signature of the a LSP. 
However, since the Mp measurement at colliders |120j . 
which would have been a decisive test of the G LSP, 
seems cosmologically disfavored, it will be a challenge 
to distinguish between the a LSP and the G LSP. 

For m fl = 100 GeV and m s = 110 GeV, for ex- 
ample, the T\ lifetime in the a LSP scenario (f2TJ|) can 
range from 0(0.01 s) for f a = 5 x 10 9 GeV to O(10 h) 
for f a = 5 x 10 12 GeV. In the G LSP case, the cor- 
responding lifetime (|10p can vary over an even wider 
range, e.g., from 6 x 10~ 8 s for = 1 keV to 15 years 
for = 50 GeV. Thus, both a very short lifetime, 
t^j < ms, and a very long one, > days, will point 
to the G LSP. On the other hand, if the LSP mass can- 
not be measured kinematically and if = 0(0.01 s)— 
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0(10 h), the stau lifetime alone will not allow us to 
distinguish between the a LSP and the G LSP. 

The situation is considerably improved when one 
considers the 3-body decays ri — > rja/G. From the 
corresponding differential rates [124] . one obtains the 
differential distributions of the visible decay products. 
These are illustrated in Fig. [7] (from [124] ) in terms of 

d 2 r(r 1 -> rji) 



1 

T7i: 



dxjdcos 6 



(21) 



where x 7 = 2E J /mf 1 is the scaled photon energy, 9 
the opening angle between the directions of 7 and r, 



r(f 1 ^r 7? ;x™ t ,^ ut )^ 



l-Ai 



d 2 r(n 



dcos( 



r 7 ^j 



dx 7 d cos 



(22) 



the respective integrated 3-body decay rate with the 
cuts x 1 > a£ ut and cos# < 1 — xg ut , and = mf/m^. 
Note that (f2Tj) is independent of the 2-body decay, the 
total NLSP decay rate, and the Peccei-Quinn/Planck 
scale. The figure shows (|2ip for the axino LSP (i = a) 
with m~/m~ <C 1 (upper panel) and the gravitino 

10 MeV (lower panel), where 



LSP (t = G) with 
m fl = 100 GeV, j 



G 



110 GeV, and x c 7 ut 



0.1. In the G LSP case, the events are peaked only in 
the region where photons are soft and emitted with a 
small opening angle with respect to the tau (8 ~ 0) . In 
contrast, in the a LSP case, the events are also peaked 
in the region where the photon energy is large and the 
photon and the tau are emitted back-to-back (9 ~ n). 
Thus, if the observed number of events peaks in both 
regions, this can be evidence against the gravitino LSP 
and a hint towards the axino LSP [124] Fi 

To be specific, with 10 4 analyzed stau NLSP de- 
cays, we expect about 165±13 (stat.) events for the 
a LSP and about 100±10 (stat.) events for the G 
LSP |124j . which will be distributed over the corre- 
sponding (x 7 , cos #)-planes shown in Fig. [7] In partic- 
ular, in the region of x 1 > 0.8 and cos# < —0.3, we 
expect about 28% of the 165±13 (stat.) events in the a 
LSP case and about 1% of the 100±10 (stat.) events in 
the G LSP case. These numbers illustrate that O(10 4 ) 
of analyzed stau NLSP decays could be sufficient for 
the distinction based on the differential distributions. 
To establish the feasibility of this distinction, dedi- 
cated studies including details of the detectors and the 
additional massive material will be crucial |118j . 

4.3 Probing the Peccei-Quinn Scale f a and ma 

If a is the LSP and T\ the NLSP, the analysis of the 2- 
body decay f\ — > to, will allow us to probe the Peccei- 
Quinn scale f a and the axino mass ma. In fact, the 



18 There is a caveat: If < ma < mf l and r{j\ — > 
aX) 3> r(ri — > G X), one would still find the distribution 
shown in the upper panel of Fig. [7] The axino would then 
eventually decay into the gravitino LSP and the axion. 



1> 
to 
o 
u 




0.4 0.6 



OS 

(0 

o 
u 




0.2 0.4 0.6 0.8 1 

Fig. 7. The normalized differential distributions (|2ip of 
the visible decay products in the decays fi — > r + 7 + a/G 
for the cases of the a LSP (upper panel) and the G LSP 
(lower panel) for m ?1 = 100 GeV, \i - B, m s = 110 GeV, 
m§/m§ <C 1, and mg = 10 MeV. The cut parameters are 
set to a; 7 ut = Xg ut =0.1. The contour lines represent the 
values 0.2, 0.4, 0.6, 0.8, and 1.0, where the darker shading 
implies a higher number of events. From [124] . 



measurement of t?, (l20l) with methods described in 



Sect. 13.31 leads to the following estimate of the Peccei- 
Quinn scale f a |124j : 



2 n 2 

O. 



aYYl 10 n GeV) 2 (l 

' ™ti \ ( m B 

.100 GeV/ V 100 GeV 



V25 sJ 



(23) 



which can be confronted with f a limits from astrophys- 
ical axion studies and axion searches in the labora- 
tory pi37U1381ll391fT53] . Indeed, we expect that m fl 
and m,g will already be known from other processes 
when the t\ NLSP decays are analyzed; cf. Sect. 13.31 
The dependence on ma is negligible for ma/m^ ^0.1. 
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For larger values of ma, the f\ NLSP decays can be 
used to determine ma from the kinematics of the 2- 
body decay, i.e., from a measurement of the energy of 
the emitted tau E T , 



ma = J m 2 fi +m%- 2m fl E T , (24) 

with an error governed by the experimental uncertain- 
ties on mfj and E T . As is evident from (fl7|) and (fT8|) . 
the determination of both the Peccei-Quinn scale f a 
and the axino mass ma is crucial for insights into the 
cosmological relevance of the axino LSP. 



5 Conclusion 

Dark matter is strong evidence for physics beyond the 
Standard Model. Extending the Standard Model with 
SUSY, an electrically neutral and color neutral LSP 
becomes a dark matter candidate for conserved PL- 
parity. I have shown that the neutralino Xii the grav- 
itino G, and the axino a can be the LSP and as such 
explain the non-baryonic dark matter in our Universe. 
The neutralino Xi is already part of the MSSM which 
provides a solution of the hierarchy problem and allows 
for gauge coupling unification. Being the superpartner 
of the graviton and the gauge field associated with su- 
pergravity, the gravitino G is equally well motivated 
with a mass mg that reflects the SUSY breaking scale. 
As the superpartner of the axion, also the axino a ap- 
pears naturally once the strong CP problem is solved 
with the Peccei-Quinn mechanism in a SUSY setting. 

While mass values and interactions can be very dif- 
ferent for the Xi> C> an( i &i I have illustrated for each 
of these LSP candidates that natural regions in the 
parameter space exist in which i?LSP = fldm- These 
regions are limited by bounds from electroweak pre- 
cision observables, B-physics observables, Higgs and 
sparticle searches at LEP, and by BBN constraints. 
The constraints from Sldm and BBN also imply seri- 
ous upper limits on the reheating temperature after 
inflation Tr which can be relevant for models of infla- 
tion and baryogenesis. 

Most promising are the experimental prospects in 
the case of the x° LSP. Being a WIMP, the x°i LSP 
should be accessible in direct and indirect dark mat- 
ter searches. Indeed, first hints might have already 
been found in the EGRET data [MEI]. With ongo- 
ing indirect searches, the increasing sensitivity of di- 
rect searches, and the advent of the LHC at which 
Xi dark matter could be produced, we will be able to 
test whether these hints are indeed the first evidence 
for the existence of SUSY dark matter. While an ex- 
cess in missing transverse energy is expected to be the 
first evidence for SUSY at the LHC already within 
the next three years, the identification of the x? being 
the LSP will require the reconstruction of the SUSY 
model realized in nature. If superparticles are within 
the kinematical reach, precision studies at the ILC will 
be crucial for this endeavor. 



In the G/a LSP scenarios with conserved R-parity, 
no dark matter signal should appear in direct or in- 
direct searches. However, since an electrically charged 
LOSP such as the t\ is viable in the G/a LSP scenar- 
ios, (quasi-) stable ti's might occur as muon-like par- 
ticles instead of an excess in missing transverse energy. 
Indeed, an excess of (quasi-) stable ri's could appear 
as an alternative first evidence for SUSY at the LHC in 
the next three years. Because of the severe upper limits 
on the abundance of stable charged particles [4], one 
would then expect that the f\ is the NLSP that decays 
eventually into the G/a LSP or that R-parity is bro- 
ken. A distinction between these scenarios will require 
the analysis of the t\ decays. For this challenge, the 
ILC with its tunable beam energy seems crucial [1 151 
HTBIIT5T1ITT81IT551 . 

Table [2] presents an overview of the SUSY dark 
matter candidates discussed in this review. As the LSP, 
each of them — the lightest neutralino x? > the gravitino 
G, or the axino a — could provide fidm and could be 
produced and identified at colliders in the near future. 

I would like to thank the organizers of SUSY 2007 
for inviting me to an exciting and stimulating con- 
ference. I am grateful to A. Brandenburg, L. Covi, 
A. Freitas, K. Hamaguchi, G. Panotopoulos, T. Plehn, 
J. Pradler, L. Roszkowski, S. Schilling, N. Tajuddin, 
Y. Y. Y. Wong, D. Wyler, and M. Zagermann for valu- 
able discussions and collaborations on the topics cov- 
ered in this review. 
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